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ABSTRACT 



A quantitative study of the observable radio signatures of the sausage, kink, and 
torsional MHD oscillation modes in flaring coronal loops is performed. Considering 
first non-zero order effect of these various MHD oscillation modes on the radio source 
parameters such as magnetic field, line of sight, plasma density and temperature, elec- 
tron distribution function, and the source dimensions, we compute time dependent radio 
emission (spectra and light curves). The radio light curves (of both flux density and de- 
gree of polarization) at all considered radio frequencies are than quantified in both time 
domain (via computation of the full modulation amplitude as a function of frequency) 
and in Fourier domain (oscillation spectra, phases, and partial modulation amplitude) 
to form the signatures specific to a particular oscillation mode and /or source parameter 
regime. We found that the parameter regime and the involved MHD mode can indeed 
be distinguished using the quantitative measures derived in the modeling. We apply 
the developed approach to analyze radio burst recorded by Owens Valley Solar Array 
and report possible detection of the sausage mode oscillation in one (partly occulted) 
flare and kink or torsional oscillations in another flare. 

Subject headings: Sun: flares — Sun: oscillations — Sun: radio radiation — Sun: X-rays, 
gamma rays — Sun: corona — (magnetohydrodynamics:) MHD 



1. Introduction 



Oscillations and quasi- periodic pulsations (QPPs) in radio emission from solar flares have been 



observed for half a century (lYoung et al.lll96ll ). Pulsation phenomena can be formally categorized 
by wavelength of the emission and by frequency /period of the QPPs. For example, at meter 
and decimeter wavelengths, rapid QPPs occur with periods of tens to hundreds of milliseconds. 
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Chen et al.l 120111 ). Longer-period QPPs, with periods r > 10 s are observed in microwaves : 



such QPPs are ofte r i accompanied by QPPs in hard X-rays (HXR; e.g.. iParks &: Wincklei 
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Numerous models have been proposed to account for th e various types of oscillations, QPPs, 



and r elated phenomeiia obs erved at radio wavelengths (e.g., lAschwandenl 119871 : iNakariakov et al 
20031 : [Fleishman et al.l 120081 . and references therein) . Generally speaking, for radio emissions pro- 
duced by coherent radiation processes (such as fast millisecond pulsations), QPPs are believed to 
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Korsakov &: Fl eish manlll998l ). In contrast, for radio emission produced by incoherent gyrosyn 



chrotron (GS) radiation from energetic electrons, QPPs are believed to result from modulation of 
the source parameters (such as the energetic electron distribution, the magnetic field strength, the 
line of sight, etc.) v ia magnetohydrodynamic (MHD) oscillations (kink, sausage, or torsional modes) 
(jAschwandenll2004l ) and/or modulation of electron distribution function via time-depe ndent accel- 



eration and injection. Most of the reported so far Q PPs of GS microwave emission (e.g.. lKane et al 



1983 



Nakaiima et al.lll983l : [Fleishman et al.ll2008l ) are consistent with modulation of fast electron 



acceleration/injection, while there is no compelling case in which modulat ion of GS emission by a n 
MHD oscillation mode would be firmly confirmed. As has been shown by lFleishman et al.l (|2008l ). 
a well spectrally resolved radio spectrum is crucial to distinguish between MHD oscillations and 
modulation of fast electron injection/acceleration. Another component needed to firmly detect a 
particular MHD oscillation mode from the observed radio data is detailed quantitative modeling 
of observable measures, such as oscillation phase or modulation power, of the GS flaring emission 
modulated by one or another MHD mode, which is yet unavailable. 

This paper attempts to partly remedy the situation by analyzing dependence of GS radia- 
tion on modulation of the source parameters by various MHD modes, namely, large-scale kink, 
sausage, and torsional oscillation modes, while we do not consider a special case of the longitudinal 
slow mode (essentially, the sound wave) where only the plasma density oscillates since this mode 
can only affect low-freq u ency radiation from a dense plasma, which has already been studied by 
Nakariakov &: Melnikovl (j2006l ). For each mode, time-domain and Fourier analyses are performed 
on the radio emission and polarization to characterize the properties of the pulsations, and the 
results are evaluated to identity a signature by which each oscillation mode can be distinguished 
in observational data. In contrast to approximate sausage-mode modulation modeling performed 
by [Fleishman et al.l (|2008l ) based on relatively crude Dulk-Marsh approximation, in this study the 
mi crowave (GS plus free-free) ein ission is modeled using precise fast GS codes recently developed 
by lFleishman fc Kuznetsovi (|2010l ). 



The results of the modeling are then used to search for the corresponding QPP signatures in the 
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total power microwave database available from the Owens Valley Solar Array (OVSA). Specifically, 
we looked through well-calibr ated data record ed during 2001-2002 when microwave emission from 
412 solar flares was recorded (|Nita et al.ll2004l ). We found the the vast majority of the radio bursts 
does not display unambiguous QPPs in the dynamic spectra; only about 10% of all events show 
time variations at the light curves that can be classified as QPPs. From them we selected a few 
events in which the presence of QPPs is apparent from their dynamic spectra; for those selected 
events we computed all possible quantitative measures to be compared with the corresponding 
measures derived from the modeling in the attempt to detect an indication of one or another MHD 
loop oscillation mode. In most of the case we were unable to firmly identify the genuine cause of 
the QPPs; however, in two best cases, presented in the paper, the observed QPPs indeed could 
be produced by an MHD mode: the sausage mode in an (occulted) flare of 28 August 2002 21:41 
UT and kink or torsional mode in a flare of 20 July 2002 21:27 UT. We note that some MHD 
oscillation modes may or may not produce corresponding oscillations in HXR and/or soft X-rays 
(SXR), which we also tested for the mentioned events using HXR data available from the Reuven 
Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and SXR data from the Geostationary 
Operational Environmental Satellites (GOES). 



2. Modeling of GS emission modulated by MHD modes 

2.1. MHD loop oscillation modes 

The different modes of MHD loop oscillations are modeled by oscillating perturbations of mag- 
netic field vector and further calculating the linear effect of these oscillations on relevant parameters 
affecting the GS radiation. In particular, the magnetic flux conservation affects cross-sectional area 
and volume of the loop. Then, assumed particle number conservation has an effect on the number 
density. The first adiabatic invariant results in oscillations of the particle energy. Oscillations of 
the magnetic vector direction implies oscillations of the viewing angle between the magnetic field 
lines and the line of sight. 

For the sake of simpli ci ty, the flariii g coro nal loop is approximated to be a cylinder as in 



Zaitsev Sz Stepanovl (jl982bl ) ; lAschwandenl (j2006l ) . The oscillation modes of a finite cylinder must 



necessarily display both time and spatial dependence; however, the latter strongly complicates the 
entire treatment as it implies building a nonun iform 3D model and so lving the radiation transfer 
equation along many lines of sight (see, e.g., Kuznetsov et al.l boill ). which is far beyond our 



intention to address a spatially uniform problem. Thus, for the sake of simplicity, we entirely 
discard any dependence of the oscillations on the spatial coordinates s, r, and if and retain only 
time dependence of the MHD modes; so for the case of the time-oscillating cross-sectional area and 
surface area, both the side and base of the cylinder were considered as affecting GS radiation by 
taking the visible area to be 
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A{t) ^ S{t) sm{e{t)) + -L\t) cos{9{t)) (1) 

where S{t) is the cyhnder side area, L{t) is the cyhnder base diameter, and 9{t) is the viewing angle 
between the cyhnder axes and the hne of sight. Accordingly, the source depth d is the smallest of 
S{t)/L{t) and 

d{t) ^ L{t)/sm{e{t)). (2) 

The large-scale (global) sausage mode is modelled as an oscillation of the magnetic field parallel 
to the coronal loop: 

B,{t) = Bo{l + asm{ujt + <p)) . (3) 

The first adiabatic invariant implies corresponding oscillations of the transverse particle momentum 
p±. Without energy losses and angular scattering of the particles, these p± oscillations imply 
variations of both particle energy and pitch-angle, i.e., variation of electron energy and angular 
distributions. If, however, electron isotropization due to angular scattering occurs faster than 
the oscillation period then the transverse energy will be redistributed between the transverse and 
parallel motions to retain the angular distribution isotropic independently on the oscillation phase 
still conserving the total energy transmitted from the oscillation to the particles. This energy 
conservation can, however, break down if the processes of the energy loss/gain are important. For 
example, the energy can decrease due to thermal conduction or radiative losses or increase due to 
a flare energy release including the Alfven pump effect driven by the very same MHD oscillations. 
Below we adopt that a balance between the energy gain and loss results in a gradual variation 
of the plasma mean (over the MHD period) temperature on top of which the MHD-wave induced 
oscillations are superimposed; thus, the temperature oscillates as 

around a mean level, T(0), w hich, perhaps, itself smoothly varies with time. The upper and lower 



energy cutoffs (-Emim -^max) (IFleishman &: Kuznetsovll2010l ) for accelerated electrons fluctuate by 



the same law. By flux conservation, this induces an oscillation of the cylinder cross-sectional area 
directly proportional to {B/Bq)~^ and, therefore, the cylinder radius ~ {B / Bq)~^/'^ , and then 
in number density of both thermal and nonthermal particles ~ B/Bq. The viewing angle is not 
affected in the sausage mode. Thus, the source visible area oscillates according to Eqn ([1]) due to 
oscillations of the cylinder dimensions, while the viewing angle remains constant and does affect 
the oscillations. 

In the kink mode, the oscillation of the magnetic field was taken as 

B^{t) = B^osm{ujt + (f)) (5) 



where is the magnetic field in the direction of the kink perturbation. The kink mode is known 
to be weakly compressible, to the order of r/A, where A is the oscillation wavelength, so there 
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are some variations of the density and the cross-sectional area in this wave mode. However, the 
adopted spatially uniform oscillations imply r/A <^ 1, which allows discarding these variations, so 
we assume a constant cross-sectional area of the cylinder in the kink oscillations; however, visible 
source area ([T|) and its depth ([2]) are affected by the magnetic field oscillations described by Eqn ([5]) 
due to the viewing angle variations. Thus, the main linear effect in these oscillations is produced 
by the viewing angle, while the magnetic field value experience only a second order oscillations. If 
the unit vector toward the observer is / = then 



IB 



cos 6{t) = = cos 00 + a sm{ujt + (j)) sin 9q cos ip, 

IB' 



(6) 



where a = B±q/Bq and ip is the angle between the direction of oscillation of B± and the line of 
sight. Therefore, there are two main extreme cases to be considered here, cos(^ = 1 and cos 99 = 0. 
In the first case, the radiation variation arises mostly from the oscillation of the viewing angle. In 
the second case, however, the viewing angle will not oscillate, and the radiation oscillation will arise 
from a second-order oscillation of the magnetic field strength. Apparently, all intermediate cases, 
< cos ip < 1 are also possible; the results will be between the mentioned extreme cases. 

Similarly, the torsional mode only affects the pitch angle of the magnetic field in the linear 
approximation. The magnetic field oscillates as 



B^{t) = B^o sm{ujt + 



(7) 



so only the viewing angle oscillates as (|Tapping||l983l ) 



cos 9{t) = cos 00 cos — sin sin (V'(i)) , 



(8) 



where ^{t) = asm cot is the pitch-angle of the magnetic field. In contrast to the kink mode, there 
no oscillations of the source depth and visible area in the torsional mode, because oscillations of 
the azimuth component B^p do not imply spatial variation of the cylinder as a whole. 



These oscillatio n modes have been studied and described ex t ensively bv iRosenbera (119701) 



Roberts et al.l (|l983l ): lAschwandenl (119871 ) : lAschwanden et al.l (|l999l ): lNakariakov et all (120031 ) 
(|2003l ). The kink and sausage modes both have fast and slow MHD wave versions. For the fast 
sausage mode, the observed period is defined by the cylinder length / and the wave phase speed Vph 



Grechnev et al, 



Tsausage ~ < 120 ^- SeCOnds 

Vph B 



(9) 



where rg is the tube radius in 10^ cm, nio is the tube density in 10^" cm ^, and B the magnetic 
field in gauss. The fast kink mode will have a period of 



■^kink ^ 100 — - — seconds 

jB 



(10) 
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where j — 1 is the number of nodes and Lg is the length of the loop in 10^ cm. The slow wave 
versions both have the same period, 

Tsiow ~ 13-^seconds (11) 



where Tq is the temperature measured in 10^ K. The sausage mode also has a n unstable propagatin g 
version with the period that is 0.6 times less fast than the fast sausage mode ([Grechnev et al.ll2003l ). 
while the torsional mode has a period 

torsional ~ IQQ V^ geCOnds, (12) 

mn 

where m is the number of windings along the loop. 



2.2. Characterization of radio oscillations 



The time-oscillating radio flux and polarization modulated by a given MHD oscillation mode 
was straightforwardly simulated using the fast GS codes: in each instance we computed an array 
of the source parameters specific to this particular mode, see § 12.11 and sent it to GS code, which 
returned the instant radio spectrum and polarization at the selected spectral range, typically, from 
1 to 100 GHz. Repeating this procedure many times during the oscillation period we, thus, formed 
the oscillating radio emission representing the subject of further analysis of the QPPs. 

As has beeri alrea dy said, QPPs can be characterized in both time and Fourier domains 
(jFleishman et al.ll2008l ). In the time domain, the modulation amplitude m{f) for each radio fre- 
quency can be calculated as 

T \ 1/2 

(13) 



m{f) 



where 



S\f,t)dt 
F{f,t)-{F{f,t)) 



S(f t)= ^■''"^'^^ '-J'''" (14) 

is the normalized modulation of the signal, / is t he frequency in GHz, t is the time in seconds, and 
F{f,t) is the total flux density at frequency / (jFleishman et al.l 120081 ). This value is a measure 
of the variation of emission with respect to the running average. The same measure can also be 
calculated for the absolute value of the degree of polarization of the emission. 

To find the oscillation periods of the modeled radio light curves, we perform a FFT of the 
GS data across all the emission frequencies. To smooth out the result of the F FT, an array of 
zero-valued points was added to the end of the data (sf.. [Fleishman et al.ll2002bl ). Since only one 
oscillation period was simulated, the Fourier coefficients of the MHD oscillations form a narrow 
peak in the FFT spectrum describing the simulated oscillation period. The radio light curves, 
however, represent a complicated nonlinear response on the MHD oscillations and so, their Fourier 
spectra can strongly deviate from the original single narrow peak. 
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At each significant Fourier peak, the phase of the pulsations relative to the oscillations at a 
reference frequency /re/j that is selected somewhere in the optically thin region, is calculated as 

where 



*'^'"'' = Vto„/..)M^(/.WI 

and y) is the coefficient of the complex FFT array at emission frequency / corresponding to 
a pulsation frequency and F*{f, u) is the complex conjugate of this coefficient. 

In addition to the above, the Fourier analysis allows for consideration of the partial modulation 
of the signal by each pulsating component at a Fourier spectrum. Analysis of partial modulation is 
necessary when there is more than one significant Fourier peak in the power spectrum. Furthermore, 
different peaks may be the results of different physical processes, or of different modes of MHD 
oscillations, and a quantitative analysis of partial modulation can help discern these different causes. 
The partial modulation amplitude mp{f ) of radiation from a limite d region of Fourier har monics 



from rii to n2 covering each main Fourier peak is defined as follows (jFleishman et al.ll2008l ): 



'mpif) 



2^ E \sif,Hnm 



n=ni 



1/2 



(17) 



where Ni is the total number of the measurements in the time domain and Ntot is the total number 
of the oscillation frequencies i^(n) where the Fourier amplitudes are determined, and the factor of 
2 arises from the fact that each partial oscillation consists of two Fourier peaks, one at positive 
frequency v and one at the corresponding —v. It is worthwhile to note that adding up all the 
available Fourier components from n = to n = Ntot — 1 yields the total modulation amplitude, 
which is equivalent to that given by Eqn (jl3p because of familiar Parceval's identity. 



Microwave emission from solar fiares differs depending on the ambient plasma density, magnetic 
field, fast electron energy or angular distributions. Below we consider a few different parameter 
combinations covering main distinct regimes of the GS emission in fiares. For each considered 
parameter regime we simulate the time domain response of the corresponding radio source on all 
mentioned MHD loop oscillations as described above. Specifically, to characterize radio response on 
the MHD oscillations and determine corresponding GS emission signatures, the following measures 
were computed on both the fiux density and polarization: 

1. the full modulation amplitude as a function of emission frequency /; 

2. the Fourier spectra of the emission at each emission frequency /; 

3. the relative phase of oscillations at Fourier peaks in the power spectrum as a function of 
emission frequency /; 
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4. the partial modulation amplitudes as a function of emission frequency / for each significant 
peak in the Fourier spectrum. 



3. Modeling results 

We considered a number of distinct parameter regimes forming the microwave emission includ- 
ing the cases of the spectral peak formation by either optical thickness effect (low-density case) 
or Razin-effect (high-density case); then, we also analyzed the case of strong magnetic field when 
distinct gyroharmonics appear at the low-frequency part (1 — 10 GHz) of the microwave spectrum. 
Furthermore, we considered different angular distributions of the radiating fast electrons, namely, 
the isotropic and anisotropic loss-cone distributions; note, that the loss-cone boundary itself can 
oscillate in response to the MHD modulation. Finally, we distinguished cases when the total power 
data are analyzed (e.g., oscillations of the source area contribute to the overall modulation) and 
when we deal with a spatially resolved data from a fixed single pixel (whose area is by definition 
constant). The implied big number of the corresponding models, see overview of them in Figured] 
results apparently in a large number of illustrating plots; most of which are collected in the on-line 
supplement to the paper. In contrast, in the main part of the paper we consider in more detail 
only a lim ited number of the cas es needed for direct comparison with either available simplified 



modeling (iFleishman et al.l l2008l ) or with observational data analyzed below in the paper. We 
adopt the following default parameters: the cylinder side area Sq = 2.5 x 10^® cm^, the cylinder 
depth Lq = 2.5 x 10® cm, the background temperature Tq = 2 x 10^ K, the thermal number density 
no = 10^ cm the nonthermal number density = 0.02 x no, the magnetic field -Bo = 50 G, 
the viewing angle Oq = 45°, the electron energy spectral index 5 = 3.5 for 0.1 < < 10 MeV, the 
angular distribution is isotropic. 

We start from the sausage mode oscillations of a radio source with default-value parameters 
where the radio spectr um peak, see Figure [It, is formed by the optical thickness effect like in 



Fleishman et al.l (|2008l ). while the Razin-effect is unimportant. Figure [2] displays the oscillating 
components of the radio emission at three characteristic frequencies selected in the optically thick, 
spectrum peak, and optically thin regions of the spectrum, respectively. One can easily notice 
that although the modulation resembles a sinusoidal signal at both thick and thin frequencies, it 
departs from the sinusoidal greatly around the spectrum peak. Figure [2)3, which means that even 
small-amplitude sinusoidal MHD oscillations can result in a nonlinear response producing higher 
harmonics in the radio light curves. This conclusion is explicitly confirmed by the Fourier spectra 
of these same oscillating components given in Figure [3t.-c. Indeed, the first Fourier harmonics 
dominates the power spectra at thick (1.5 GHz) and thin (13.8 GHz) frequency, while the second 
harmonics dominates the power spectrum around the radio spectrum peak (2.2 GHz). The bottom 
raw of the same Figure displays results of the Fourier spectra analysis: the relative phase, partial, 
and full modulation amplitudes as functions of the radio frequency. 



Apparently, the results are consistent with those of the simplified modeling reported by 



- 9 - 



Fleishman et al.l (|2008l ): the oscillation phase changes by vr around the peak frequency, while the 
modulation amplitude displays a curve with a minimum around the frequency of the spectrum 
peak and with modulation at the optically thin part exceeding by a few times the modulation at 
the optically thick part . Thus, the performed here more advanced modeling confirms conclusion of 



Fleishman et al.l (|2008l ) that the QPPs analyzed by them cannot be caused by an MHD oscillation 
mode. Then, Figures [4] and [5] display the same parameters but computed for the absolute value of 
the degree of polarization, which can be used for the QPP analysis if the polarization measurements 
are available. 

The partial modulation amplitude of the degree of polarization comes to a peak at the frequency 
range where the degree of polarization changes the sign and therefore it is around zero, while 
decrease strongly at low and high frequencies. At low and high frequencies respectively, the flux is 
in phase with the degree of polarization. Figures [2] and [H 

Another kind of observations, which can help in identifying the cause of QPPs, is the imaging 
spectroscopy data taken with spatial resolution. In particular, the detailed shapes of the modulation 
amplitudes for the total power data and single-pixel view are different from each other, see on-line 
supplement, so their comparison, when available, will yield an estimate of the source visible area 
pulsations, which is highly valuable in identifying the oscillation mode. The shape of the modulation 
amplitude curve is also different for the loss-cone distribution of the fast electrons compared with the 
isotropic case, see on-line supplement; however, the difference is not that big to allow unambiguous 
distinction between the isotropic and anisotropic electron distributions at the current level of the 
available radio data. 

The presence of low-frequency peaks in the emission is an indicator of high magnetic field. In 
emission from a pulsating source, these peaks will oscillate as well. Figured^, b. As there is a phase 
shift associated with each local peak frequency and also amplitude of the base of each such peak, 
the Fourier spectra in the region of the low-frequency peaks do not exhibit consistent Fourier peaks, 
see on-line supplement, which makes both the partial modulation amplitude and pulsation phase 
fluctuating at those low frequencies. The phase repeatedly changes by vr around the gyroharmonics 
peaks. The full modulation amplitude naturally exhibits peaks at each low-frequency emission 
peak, while for high frequencies behaves similarly to the case of the default parameters; these 
low-frequency behavior makes the case of the strong magnetic field easily distinguishable from 
the default case without gyroharmonics contribution. The degree of polarization of the emission 
behaves in a similar way, see on-line supplement for more detail. 

The sausage mode modulation of a dense radio source, no = 10^'' cm~^, where the Razin effect 
is important, is clearly distinguishable from the case without the Razin effect even though the effect 
of the optical depth can still be importari t just around the spectral peak as in Figured^, where the 



Razin cutoff frequency (e.g., ISimonlll969l ) 



B smt) 

is marked by a vertical dashed line slightly below the spectrum peak frequency. The relative phase 



- 10 - 



of the flux displays no phase shift across the entire radio spectrum, i.e., oscihation of both low- 
and high- frequency hght curves occur strictly in phase. The modulation amplitude differs from 
the case of the default parameters; although it still rises at the high frequencies it does not display 
any flat region at the low frequencies, see Figure [H 

The kink mode, which we consider here only for the line-of-sight oscillations, cos(p = 1, in a 
dense plasma volume (see on-line supplement for all other parameter regimes), is easily distinguish- 
able from the sausage mode. Indeed, unlike the sausage mode, here the flux and polarization at 
a given frequency are vr out of phase with each other. The relative phase of the flux light curves 
is constant across the entire spectrum like for the sausage mode oscillations in a dense source. 
However, dependence of the modulation amplitude on frequency is strongly different from that for 
the sausage mode. Now the modulation amplitude is large at low frequencies, then decreases and 
remains relatively constant at high frequencies. 

The behavior of the kink mode oscillations is particularly interesting in case of the default 
parameters, see on-line supplement. The explicit effect of the viewing angle oscillations and cor- 
responding source visible area oscillations are counter directed at the optically thick part of the 
spectrum, which has two observable consequences: (i) there is no phase shift across the entire spec- 
trum (unlike the sausage mode, default parameters) and (ii) the modulation amplitude decreases 
toward low frequencies, while for the sausage mode case the modulation amplitude reaches a mini- 
mum around the radio flux spectrum peak. Remarkably, that with spatially resolved observations, 
the flux density from a flxed single pixel is unaffected by the visible area oscillations; as a result 
the kink mode oscillations observed from a single pixel do produce a vr phase shift in both flux and 
polarization in the region of the spectrum peak; the shape of the modulation amplitude function 
is also significantly different from that from the total power data, see on-line supplement. Thus, 
comparison of these measures for the total power and spatially resolved data is highly informative 
to distinguish the kink mode oscillations of the flaring loop. Note, that the oscillations of the source 
visible area has no effect on the optically thin part of the spectrum; in particular, it is the case for 
the entire spectrum when the Razin effect forms the spectrum peak, see above. 

The torsional mode, like the kink mode, is to the flrst order an oscillation of the viewing 
angle, although, unlike the kink mode, the source visible area and depth are not affected by the 
oscillations. For this reason, the torsional mode oscillations are distinct from the kink mode in case 
of the default parameters, while it is undistinguishable from the line-of-sight kink-mode oscillations 
at the optically thin part of the radio spectrum, in particular, in the Razin effect regime; see on-line 
supplement. 

4. Search for MHD oscillations in flaring microw^ave bursts 

The most detailed database of the solar microwave bursts recorded with reasonably high spec- 
tral resolution (40 frequencies across the range of 1 — 18 GHz) was obtained by OVSA during solar 
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cycle # 23. iNita et al.l (|2004l ) presented a thorough statistical analysis of a subset of the entire 
database composed of 412 radio bursts recorded during 2001-2002. Given the subset is sufficiently 
large and the corresponding total power data are well calibrated, we restrict our search of the 
MHD oscillation signatures to this subset; we note, however, that a number of radio bursts with 
QPPs suitable fo r a ca se study analysis was recorded after 2002, see example studied in detail by 
Fleishman et al.l hood ). 



To identify the events candidates for further quantitative analysis we manually search through- 
out the dynamic spectra and light curves of all 412 bursts and picked up those displaying QPPs. 
Surprisingly, the total number of the events-candidates is very small: less than 10 events (i.e., less 
than 2.5%) demonstrate apparent QPPs visible in the dynamic spectrum; about 40 events (~ 10%) 
display light curves with QPPs at a few frequencies. 

We performed analysis of the QPPs in those events demonstrating pulsations in the dynamic 
spectra because this means stronger signal modulation, which is apparently more suitable for the 
quantitative study, while retain all other (potentially pulsating) events for further analysis. How- 
ever, in most of the events selected for our study we were unable to firmly determine the cause of 
the observed QPPs. Below we describe two events in which detection of MHD oscillation signatures 
is likely, August 28, 2002 21:41 UT and July 20, 2002 21:27 UT. 



4.1. August 28, 2002 

The event of August 28, 2002 21:41 UT, Figure El is particula rly interesting because it was 



reported among other partly occulted events (jKrucker &: Linll2008l ). The fact that a part of the 



flaring loop is occulted is helpful for further analysis as effect of the loop nonuniformity on the radio 
emission is reduced (no contribution from the occulted part is recorded). This implies that a uniform 
source model adopted for our analysis is better justified for interpretation of the occulted flares 
compared with an average fiare. Nevertheless, even for this favorable geometry the radio spectrum 
of this flare is noticeably broader than a typical spectrum from a uniform source indicating that 
the source inhomogeneity still plays some role here. 

The Fourier spectra of the oscillating component of the radio light curves display consistent 
peaks around lis throughout the entire radio spectrum. Figure [9l although their shapes are 
not strictly identical to each other. The phase of the pulsations, however, displays neither exact 
constancy nor distinct vr shift around the radio spectrum peak. Figure [TOb. On the other hand, we 
note that the scatter of individual data points, e.g., at frequencies 10 to 15 GHz, is of the same 
order of magnitude as the range of overall phase variation over all frequencies. Thus, given overall 
weakness of the oscillating component, ±(10 — 20) sfu, implying relatively large error in the phase 
determination, we conclude that the phase does not display significant variations over the spectral 
range considered. As has been shown in § [21 the phase constancy across the spectrum is indicative 
of the Razin-effect importance in forming the radio emission, which means a relatively dense fiaring 
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loop involved in the QPPs or the kink mode oscillations, which produce constant phase in both 
high- and low- density cases. 

To distinguish between possible MHD modes capable of producing the observed QPPs, we 
address the modulation amplitude as a function of emission frequency, Figure [TOb . c. We note 
that both partial and full modulation amplitudes have similar shape, which is consistent with that 
produced by the sausage mode oscillations but deviates from those produced by the kink (for either 
high- or low- density cases) or torsional modes. We conclude that the kink and torsional modes 
are irrelevant for the QPPs from the event studied, while the sausage mode oscillations is still a 
plausible candidate. 

To further evaluate this hypothesis we recall that sausage mode oscillations stimulate corre- 
sponding oscillations of the plasma density, temperature, and fast electron distribution function, 
which must result in QPPs in both SXR and HXR emission. To test this expectation we use the 
data available from GOES and RHESSI, respectively. Figure [TTk displays the oscillating compo- 
nents of two GOES SXR channels superimposed on a optically thin radio oscillating component; 
these light curves, although not identical, do display an overall consistent behavior, which is further 
quantified by the corresponding Fourier spectra given in Figure lllb . The Fourier spectra of the 
GOES channels do display broad peaks roughly coinciding with the radio peak around 10 s period 
(note a double peak structures in both 3.4 GHz and GOES spectra around v = 0.1 Hz). How- 
ever, the spectra do not coincide precisely with each other and the high GOES channel displays an 
additional peak around r = 17 s not present in two other light curves. 

Figures [TTb. d display oscillation amplitudes and Fourier spectra for the RHESSI data. Likewise 
the SXR, the HXR oscillations show overall similarity with the radio QPPs with the highest energy 
channel, 50-100 keV, being roughly vr out of phase with the radio light curve and lower-energy 
RHESSI channels. As for the Fourier spectra, that for the highest RHESSI channel coincides 
almost identically with the Fourier spectrum of radio oscillations at a high frequency, 10.6 GHz, 
Figure [TTb. while two other, lower-energy channels, 12-25 & 25-50 keV, display a peak at somewhat 
lower frequency, v « 0.08 — 0.085 Hz and a secondary peak at ~ 0.05 Hz. We suppose that 
oscillations of thermal plasma may contribute to this additional peak and if so, it must correspond 
to the peak in the GOES Fourier spectrum at ~ 0.06 Hz, which gives us a rough estimate of the 
QPP frequency determination accuracy of bv ~ 0.01 Hz. With this accuracy, the Fourier peaks at 
V ~ 0.09 Hz and v ~ 0.08 Hz are identical, which confirms a common nature of the radio and HXR 
pulsations, while the thermal plasma (seen via GOES and lower-energy RHESSI channels) displays 
some additional QPPs (r ~ 17 s) not seen in the microwave or HXR emission. 

We conclude that the sausage mode oscillations of a dense flaring loop could be the cause of 
the observed QPPs in the August 28, 2002 21:41 UT event, although the presented evidence is not 
fully unambiguous. In contrast, a model of QPPs produced by pulsating acceleration/injection of 
the fast electrons is less consistent with the data because this model does not predict QPPs of the 
SHR emission. We note that there is no radio polarization data for this event, which, if present. 
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might further help in identifying the relevant oscillation mode. 



4.2. July 20, 2002 

The event of July 20, 2002 21:27 UT, Figure US is different from the event of August 28, 2002 
in many respects. It is a usual, not occulted, flare with relatively large radio flux up to 10^ sfu at the 
absolute peak implying higher accuracy of the derived measures characterizing QPPs. The Fourier 
spectra display two consistent peaks across the entire radio spectrum, one at ~ 22 s (i/ ~ 0.045 Hz) 
and the other at ~ 13 s (i/ w 0.08 Hz), Figure [131 The relative phases of the QPPs at both 
oscillations are almost constant across the radio spectrum. Figure [T^ . c, again indicative of the 
Razin-effect importance in the radio spectrum formation implying relatively high plasma density in 
the flaring source. The partial modulation amplitude retains a roughly similar concave-up shape in 
both oscillations. Figure [Hb. d, same as the full modulation amplitude Figure [T^. Perhaps, this 
suggests that the ~ 13 s oscillation is a higher harmonic in the Fourier decomposition of the main 
~ 22 s oscillation. Comparison of the observed and model modulation amplitudes does not favor 
the sausage mode oscillations in this case, while consistent with either torsional or line-of-sight kink 
mode oscillations, which are very hard to distinguish from each other with the developed approach. 

Those 'transverse' MHD oscillations do not involve pulsations of the loop cross-sectional area 
or volume and, accordingly, on the plasma density, temperature, and fast electron distribution. 
Therefore, we do not expect any SXR or HXR QPPs coherent with the radio QPPs. Indeed, no 
relevant oscillations is detected in the GOES light curves (not shown here), both of which are very 
smooth in this event. The HXR emission observed with RHESSI does display fluctuations with 
some peaks in the light curves being in phase with the radio light curves, while others not in phase or 
opposite phase. As a result, the main significant Fourier peak of HXR oscillations, being consistent 
at different HXR channels, ~ 0.1 Hz, does not coincide with any of the significant Fourier peak of 
the radio data. Figure [T5l suggesting that there are oscillations in the HXR emission, but likely from 
a different cause. We conclude that kink or torsional oscillations of a dense flaring loop is a possible 
cause of the radio QPPs observed from the event of July 20, 2002, while the data is more difficult 
to reconcile with the idea of fast electron QPPs within a simple single-loop source model adopted 
in this paper: the fast e lectro n pulsations would imply coherent in-phase radio and HXR QPPs, 



like in [Fleishman et al.l (|2008l ). which is not observed in the event considered here. We cannot 
exclude, however, a possibility of a more complex source geometry in which less straightforward 
relationships between radio and X-ray QPPs are well possible. To address this question, a detailed 
imaging spectroscopy data is needed at a currently unavailable level. 
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5. Discussion 

We have described a detailed modeling of the MHD loop oscillation effect on the GS microwave 
emission from solar flares. Unlike most of other studies concentrating on observed QPP periods 
and their possible correspondence to one or another MHD mode, we address the question how 
the oscillation phase and corresponding modulation amplitude behave as a function of the radio 
frequency, give n that these QPP me asures were shown to be highly efficient to identify the true 



cause of QPPs (jFleishman et al.ll2008l ). We found that the QPP relative phase behavior is distinctly 
different for the cases of relatively strong magnetic field when the phase repeatedly changes at each 
of the local peaks at low gyroharmonics and of either dense or tenuous flaring loops with a lower 
field, displaying no or vr phase shift around the (single) spectral peak, respectively. Furthermore, 
the sausage mode oscillations are distinguishable from the kink and torsional mode oscillations 
using the frequency behavior of the modulation amplitude. The kink and torsional modes are 
distinguishable from each other in the case of tenuous source, when oscillations of the source visible 
area dominate the kink mode effect on the radio QPPs in the optically thick region. However, the 
kink and torsional mode oscillations in a dense loop are indistinguishable from each other as far 
as linear order effect (which is identical for these modes) dominates the radio emission response on 
these MHD oscillations. 

We have used the modeling results to analyze the OVSA microwave burst database and search 
for the expected MHD oscillation signatures in some bursts displaying QPPs. Although we did 
not find many promising events-candidates, we did identified two reasonably good examples in 
which detection of one or another MHD oscillation mode is likely. Specifically, in one of them 
the radio QPPs might be caused by the sausage mode, while in the other one by either kink or 
torsional mode. Comparison of the radio QPPs with SXR and HXR QPPs favors the MHD loop 
oscillations for these two events compared with t he fast electron injectio n/acceleration oscillations 



detected in another event using a similar analysis (jFleishman et al.ll2008l ). However, this conclusion 



must be taken with some caution (given incompleteness and insufficient accuracy of the data) and 
must be further evaluated and tested on other events with radio QPPs; we plan to perform a more 
comprehensive statistical study based on full available OVSA database elsewhere including analysis 
of the polarization data when available, because as evident from the above modeling, such a study 
would benefit greatly from accurate polarization data. 

We anticipate that similar data analysis will be even more important for soon-to-be-available 
imaging spectroscopy data after on-going OVSA expansion (EOVSA) has been completed. In 
particular, all EOVSA antennas will be equipped with dual-polarization feeds for simultaneous 
right-hand and left-hand circular polarization measurements. Furthermore, the expansion will 
greatly improve the image quality and the frequency resolution, which are all needed to perform a 
more meaningful analysis and yield more reliable conclusions on the cause of radio QPPs in each 
particular case. 
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Fig. 1. — : Flux density (left column) and polarization (right column) for, from top to bottom, 
examples of dynamic spectra with the low gyroharmonics present, and spectra for default param- 
eters, single pixel emission, anisotropic pitch-angle distribution, Razin effect parameters (dashed 
line marks Razin cutoff frequency), and low-frequency spectrum peaks. 
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(a) 1.5 GHz (b) 2.2 GHz (c) 13.8 GHz 

Fig. 2. — : Sausage mode, default parameters, oscillation of flux at indicated frequencies. Note that 
(i) the low- and high- frequency light curves are out of phase and (ii) the oscillating component 
deviates strongly from a sinusoidal around the spectrum peak frequency, e.g., 2.2 GHz. 
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Fig. 3. — : Sausage mode, default parameters, oscillations of the flux density, (a)-(c) Fourier 
spectra for the emission frequencies indicated, vertical lines mark off the peak regions considered 
for partial modulation according to Eqn (|17p : note enhanced second harmonics at 2.2 GHz. (d) 
Relative phase of flux, solid line is the fundamental harmonics (1 Hz), dashed line is the second 
harmonics (2 Hz); note the vr phase change around the spectrum peak frequency for both oscillation 
harmonics, (e) Partial modulation amplitude for the indicated range of Fourier coefficients for two 
harmonics: the solid (dashed) line represents integration between the solid (dashed) vertical lines 
in panels (a)-(c). (f) Full modulation amplitude for the oscillating component; it has a deep 
minimum around the spectrum peak frequency likewise the partial modulation amplitude of the 
first harmonics; in contrast, the second harmonics modulation amplitude has a maximum at this 
range, see panel (e). 
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Fig. 4. — : Sausage mode, default parameters, oscillation of polarization at indicated frequencies. 
The oscillating components at low frequencies are vr out of phase relative to the high frequencies. 
The high- (low-) frequency polarization is in phase with the high- (low-) frequency flux density, see 
Figure M 
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Fig. 5. — : Sausage mode, default parameters, oscillations of polarization, (a)-(c) Fourier spectra 
for the emission frequencies indicated. Vertical lines mark off peaks considered for partial modula- 
tion, like in Figure [3l (d) Relative phase of flux. Solid line is the fundamental harmonic, dashed 
line is the second harmonic, (e) Partial modulation amplitude for the indicated range of Fourier 
coefficients for both harmonics, (f) Full modulation amplitude for both harmonics. The modu- 
lation amplitudes reach the peak values around the spectrum peak frequency, where the sense of 
polarization changes from X-mode (at high frequencies) to 0-mode (at low frequencies). 
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Fig. 6. — : Sausage mode, Razin effect, oscillations of flux and polarization: (a)-(f) same as in 
Figures [3] and (g)-(l) same as[5l 
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Fig. 8. — : August 28, 2002. Top: Dynamic spectrum of the flux density normalized by the absolute 
peak flux of 212 sfu. The vertical stripe-like features in the dynamic spectrum are indicative of 
prominent QPPs in the event. Vertical dashed lines indicate examined region. Bottom: oscillating 
components at frequencies 5.4 (dashed) and 10.6 GHz (solid), shown in the top panel by the dashed 
and solid lines, respectively. The low- and high- frequency oscillations occur roughly in phase. 
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Fig. 9. — : August 28, 2002. Fourier spectra for indicated frequencies. Dashed lines indicate peak 
regions considered for further partial modulation calculation, see Figure [TOl below. 
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Fig. 10.—: August 28, 2002. (a): Relative phase of flux oscillations at 0.09 Hz. (b): Partial 
modulation amplitude at the harmonic range selected for integration by the dashed lines in Figure [9j 
(c): Full modulation amplitude during the time segment selected for the analysis; see dashed 
vertical lines in Figure [H top. 



-27- 




21:43:00 21:43:20 21:43:40 

Start Time (28-Aug-02 21:42:50) 



21 :44:00 

a 



1.0 






/■\ 






0.8 






/' \ 
/ \ 
/' \ 


i 




0.6 






; 


1 

1 


/ 'A 'i - 


0.4 








\ 1 




0.2 




.-/^ /~'\ i 

y \ / \w- 




\i 




0.0 












0.00 


0.02 0.04 


0.06 


0.08 


0.10 0.12 0. 



Frequency, Hz 



0. 0-2 

X) 

^ 0.1 

Q_ 

E 

D 

^ 0.0 
_o 

.1 -0.1 




-0.2 

21:43:00 21:43:10 21:43:20 21:43:30 21:43:40 21:43:50 21:44:00 
Stort Time (28-Aug-02 21:43:00) 

C 



1 .0 
0.8 






/ 

/ 


\ 


,> - 


0.6 


- 




/ 

; 


\ 

\ 




0.4 




/ 




\ 


/ w \ \ - 








/ 


^'"^ / 


1 '\ \ - 


0.2 






X-- / 


\ X 


0.0 













0.00 0.02 0.04 0.06 0.08 0.10 0.12 
Frequency, Hz 



0.14 



Fig. 11. — : August 28, 2002. (a): Dashed and dash-dotted hues show the oscillating components 
of the low- and high- energy GOES channels, respectively; the solid line shows the radio light curve 
at 3.4 GHz for comparison, (b): the corresponding Fourier spectra for the considered oscillating 
components. Note overall similarity of the broad double peaks around 0.1 Hz in the Fourier spectra 
of both SXR channels and the radio light curve, (c): Dashed, single-dotted, and triple-dotted 
lines show oscillation components for the RHESSI channels 12-25 keV, 25-50 keV and 50-100 keV 
respectively. The solid line represents the radio emission at 10.6 GHz, which is roughly tt out of 
phase compared with the highest RHESSI channel, 50-100 keV. (d): Fourier spectra of the above 
channels. Note that the solid (radio) Fourier spectrum is almost identical to the triple-dotted 
(highest energy RHESSI channel) Fourier spectrum. 
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Fig. 12. — : July 20, 2002. Top: Dynamic spectrum of the flux density normalized by the abso- 
lute peak flux of 9858 sfu. Vertical dashed lines indicate examined region. Bottom: Oscihating 
components for two radio frequencies, 10.6 and 16.4 GHz, shown by the dashed and solid hues, 
respectively. 
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Fig. 13. — : July 20, 2002. Fourier spectra for indicated frequencies. Vertical dashed lines indicate 
the peak regions considered for calculation of partial modulation amplitudes. 
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Fig. 14. — : July 20, 2002. Top: Phase and partial modulation amplitude for the Fourier peak at 
~ 22 s (iz = 0.045 Hz). Middle: Phase and partial modulation amplitutde for the Fourier peak 
at ~ 13 s (i/ = 0.077 Hz). Bottom: Full modulation amplitude shown here in linear scale giving 
better view of the modulation amplitude variation with frequency. General similarity of the partial 
modulation amplitudes for the two Fourier peaks and of the full modulation amplitude is suggestive 
of these two oscillations to be driven by the same mechanism. 
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Fig. 15. — : July 20, 2002. Top: Dashed, single-dotted, and triple-dotted lines show oscillation 
components for the RHESSI channels 12-25 keV, 25-50 keV, and 50-100 keV respectively. The solid 
line represents the radio emission at 16.4 GHz. Bottom: Fourier spectra of the above channels. 
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On-line Album of quasiperiodic pulsations (QPP) of flaring radio emission driven by 

MHD loop oscillations 

Like in the main part of the paper, we adopt the fohowing default parameters: the cyhnder 
side area 5*0 = 2.5 x 10^® cm^, the cyhnder depth Lq = 2.5 x 10^ cm, the backgromid temperature 
To = 2 X lO'^ K, the thermal number density no = 10^ cm~^, the nonthermal number density 
nf, = 0.02 X no, the magnetic field -Bo = 50 G, the viewing angle ^o = 45°, the electron energy 
spectral index S = 3.5 for 0.1 < < 10 MeV, the angular distribution is isotropic, unless otherwise 
is explicitly stated. For the Razin-effect parameter regime the thermal number density is no = 
10^*^ cm~^. For the anisotropic distribution of the fast electrons the loss-cone angle 9 = 60°. For 
the high magnetic field case the magnetic field Bq = 500 G. 

A. Sausage mode 

In addition to the plots for the default-value parameters and Razin-effect regime presented 
in the paper itself, we give here corresponding plots for the single pixel view when the visible 
area A remains fixed, while all other parameters oscillate as before, Figure [TBI as well as the case 
of anisotropic loss-cone distribution of fast electrons, when the loss-cone boundary oscillates in 
addition to all other source parameters Figure [171 

The presence of low-frequency peaks in the emission is an indicator of high magnetic field. 
In emission from a pulsating source, these peaks oscillate as well, Figure [TSl accordingly, the 
Fourier spectra consist of many peaks, Figure [T9l The oscillation phase repeatedly changes around 
the gyroharmonics and the modulation amplitude fluctuates at the low-frequency spectral range. 
Figure [20l Figures [2TVI23I demonstrate the same parameters for the degree of polarization, while 
Figures [2411291 display the same parameters for a different viewing angle, 80°. 

B. Kink mode 

The kink mode, default parameters. Figures [30H33l is distinguishable from the sausage mode. 
Indeed, in the optically thin regime at high frequencies, the flux and polarization are vr out of phase 
with each other, see Figures [30] and [321 the oscillation phase does not change with frequency, and 
the modulation amplitude curve has a distinctly different shape decreasing at lower frequencies (at 
the optically thick region). Figure [311 

The single pixel view is significantly different from the total power data because no visible 
source area oscillations are involved in the QPPs from a single pixel. Figures [M14371 Here the phase 
experiences the vr shift around the spectrum peak frequency and the modulation amplitude forms a 
curve with a clear minimum around the spectrum peak. These differences between the total power 
and single pixel view results could be highly valuable to distinguish the kink mode in the imaging 
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spectroscopy observations. 

Figures [38H1I] display the results for the kink mode oscillations in the plane transverse to the 
line of sight. In this case the effect is second order over the magnetic field perturbation. Accordingly, 
the second harmonics oscillation dominates the Fourier spectrum and the modulation amplitudes 
are very small. 



C. Torsional mode 

The torsional mode oscillations are to first order the oscillations of the viewing angle with no 
source visible area or depth oscillations. As such, the torsional mode produces results similar to 
the kink mode line-of-sight oscillations for the single pixel view and for the Razin-effect parameter 
regime, see Figures [53ti6^ 
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Fig. 16. — : Sausage mode, default parameters, single pixel view, (a)-(c) Fourier spectra of the 
flux density oscillating components for the emission frequencies indicated. Vertical lines mark off 
peaks considered for partial modulation, (d) Relative phase of flux. Solid line is the fundamental 
harmonics (1 Hz), dashed line is the second harmonics (2 Hz), (e) Partial modulation amplitude 
for the indicated range of Fourier coefficients for two harmonics: the solid (dashed) line repre- 
sents integration between the solid (dashed) vertical lines in panels (a)-(c). (f) Pull modulation 
amplitude, (g)-(l) Same for polarization. 
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Fig. 17. — : Same as in Figure [16] for sausage mode, default parameters but with anisotropic 
loss-cone pitch-angle distribution; 9c = 60°. 
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Fig. 18. — : Sausage mode, high magnetic field, the viewing angle = 45°: low frequency gyro- 
harmonics are present in the spectrum. Top: Dynamic spectrum of simulation - oscillations of 
peaks are clearly visible. Middle: Locations of low-frequency peaks over one period of oscillation. 
Bottom: Variation of the peak positions over one period of oscillation. 
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Fig. 19. — : Sausage mode, high magnetic field. Fourier spectra at a few selected radio frequencies 
selected in the low-frequency range. Strong departure from a simple sinusoidal oscillations is clearly 
visible. 
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Fig. 20. — : Sausage mode, high magnetic field. Left: Relative phase of the flux density at = 1 Hz. 
Right: Pull modulation amplitude of the flux density. For higher frequencies, it is similar to the 
cause of default parameters. At low frequencies, however, the modulation amplitude fluctuates 
strongly: it increases noticeably at the frequencies where the gyroharmonics are enhanced. 
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Fig. 21. — : Sausage mode, high magnetic field, oscillations of polarization. Same as in Figure [THl 
but for the polarization. 



-40- 




° 0.02 



0.08 
0.06 \ 



CD 

I 0.04 



0.02 
0.00 







5 10 15 20 25 
Frequency, Hz 

(a) 1.6 GHz 




5 10 15 20 25 
Frequency, Hz 



(c) 1.9 GHz 



CD 
O 



0.001 5 



0.001 - 



0.0005 



0.0000 



o 



0.01 5 



0.010 



0.005 



0.000 




5 10 15 20 25 
Frequency, Hz 

(b) 1.7 GHz 




0) 0.08 

^ 0.06 

0.04 

0.02 
0.00 








5 10 15 20 25 
Frequency, Hz 

(d) 2.1 GHz 



5 10 15 20 25 
Frequency, Hz 

(e) 2.3 GHz 

Fig. 22. — : Sausage mode, high magnetic field, oscillations of polarization. Same as in Figure [19] 
but for the polarization. 
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Fig. 23. — : Sausage mode, high magnetic field, oscillations of polarization. Same as in Figure [20l 
For higher frequencies, the modulation amplitude is similar to the cause of default parameters. At 
low frequencies, it fluctuates decreasing between gyroharmonics emission frequencies. 
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Fig. 24. — : Sausage mode, high magnetic field, oscillations of the flux density. Same as in Figure [TSl 
but for viewing angle 80°. 
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Fig. 25. — : Sausage mode, high magnetic field, oscillations of the flux density. Same as in Figure [T9l 
but for viewing angle 80°. 
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Fig. 26. — : Sausage mode, high magnetic field, oscillations of the flux density. Same as in Figure [20l 
but for viewing angle 80°. 
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Fig. 27. — : Sausage mode, high magnetic field, oscillations of the polarization. Same as in Figure [2T] 
but for viewing angle 80°. 
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Fig. 28. — : Sausage mode, high magnetic field, oscillations of the polarization. Same as in Figure [22] 
but for viewing angle 80°. 
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Fig. 29. — : Sausage mode, high magnetic field, oscillations of the polarization. Same as in Figure[ 
but for viewing angle 80°. 
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Fig. 30. — : Kink mode, default parameters, oscillation of the flux density at indicated frequencies. 
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Fig. 31. — : Kink mode, default parameters, oscillations of the flux density, (a)-(c) Fourier spectra 
for the emission frequencies indicated, vertical lines mark off the peak regions considered for partial 
modulation, (d) Relative phase of flux, solid line is the fundamental harmonics (1 Hz), dashed 
line is the second harmonics (2 Hz); note the phase constancy across the entire spectral range, (e) 
Partial modulation amplitude for the indicated range of Fourier coefHcients for two harmonics: the 
solid (dashed) line represents integration between the solid (dashed) vertical lines in panels (a)-(c). 
(f ) Full modulation amplitude for the oscillating component; it decreases toward the optically thick 
part of the spectrum in contrast to sausage mode, which displays a minimum modulation around 
the spectrum peak frequency. 
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Fig. 32. — : Kink Mode, default parameters, oscillation of polarization at indicated frequencies. 
Note, the polarization is n out of phase with the flux density at high frequencies, while in phase at 
low (optically thick) frequencies, sf Figure [30l 
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Fig. 33. — : Kink mode, default parameters. Same as in Figure[33]but for oscillations of polarization. 
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Fig. 34. — : Kink mode, single pixel view, oscillation of flux at indicated frequencies. 
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Fig. 35. — : Kink mode, single pixel view, oscillations of flux, (a)-(f) Same as in Figure [3T1 Note 
that unlike Figure [3T| here the phase experiences the vr shift, panel (d), and modulation amplitude 
displays a curve with a deep minimum, panels (e, f), around the spectrum peak frequency 
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Fig. 36. — : Kink Mode, single pixel view, oscillation of polarization at indicated frequencies. Now 
at any frequency (thin or thick) the polarization is tt out of phase relative to the flux density at 
the same frequency. 
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Fig. 37. — : Kink mode, single pixel view, oscillations of polarization. (a)-(c)Fourier spectrum for 
the emission frequency indicated. Dashed lines mark off peaks considered for partial modulation. 

(d) Relative phase of flux. Solid line is the fundamental harmonic, dashed line is second harmonic. 

(e) Partial modulation amplitude for the indicated range of Fourier coefficients, for the first two 
harmonics, (f) Full modulation amplitude. 
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Fig. 38. — : Kink mode, no line-of-sight oscillations, oscillation of the flux density at indicated 
frequencies. 
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Fig. 39. — : Kink mode, no line-of-sight oscillations, oscillations of flux, (a)-(f) Same as in Fig- 
ure [311 Note enhanced second harmonic oscillations at all frequencies and overall small modulation 
amplitude as it is a second-order effect in this case. 
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Fig. 40. — : Kink Mode, no line-of-sight oscillations, oscillation of polarization at indicated frequen- 
cies. 
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Fig. 41. — : Kink mode, no line-of-sight oscillations, oscillations of polarization. (a)-(f)Same as in 
Figure [33l 
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Fig. 42. — : Kink mode, high magnetic field. Same as in Figure [TSl 
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Fig. 43. — : Kink mode, low-frequency peaks in the flux density. Same as in Figure [T9l 
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Fig. 44. — : Kink mode, low-frequency peaks in the flux density. Left: Relative phase of flux. Right: 
Full modulation amplitude of flux. Note oscillating behavior of the modulation amplitude at low 
frequencies. 
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Fig. 45. — : Kink mode, high magnetic field, oscihations of polarization. Same as in Figure [211 
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Fig. 46. — : Kink mode, high magnetic field, oscihations of polarization. Same as in Figure 
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Fig. 47. — : Kink mode, high magnetic field, osciUations of polarization. Same as in Figure 
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Fig. 48. — : Kink mode, high magnetic field, osciUations of the fiux density, viewing angle 80°. 
Same as in Figure [IHl 
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Fig. 49. — : Kink mode, high magnetic field, oscihations of the flux density, viewing angle 80°. 
Same as in Figure [191 
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Fig. 50. — : Kink mode, high magnetic field, oscihations of the flux density, viewing angle 80°. 
Same as in Figure [20l Note prominent oscillating behavior of the modulation amplitude at low 
frequencies. 
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Fig. 51. — : Kink mode, high magnetic field, oscillations of polarization, viewing angle 80°. Same 
as in Figure EU 
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Fig. 52. — : Kink mode, high magnetic field, oscillations of polarization, viewing angle 80°. Same 
as in Figure [22l 
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Fig. 53. — : Kink mode, high magnetic field, oscillations of polarization, viewing angle 80°. Same 
as in Figure [23l 
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Fig. 54. — : Torsional mode, default parameters, oscillation of the flux density at the indicated 
frequencies. 
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Fig. 55. — : Torsional mode, default parameters, oscillations of the flux density, (a)-(f) Same as in 
Figure [3TJ 
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Fig. 56. — : Torsional mode, default parameters, oscillation of polarization at the indicated frequen- 



cies. 
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Fig. 57.—: 
Figure [33l 



Torsional mode, default parameters, oscillations of polarization, (a)-(f) Same as in 
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Fig. 58.—: 
Figure [TBI 
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Torsional mode, Razin effect, oscillations of flux and polarization, (a)-(l) Same as in 
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Fig. 59. — : Torsional mode, high magnetic field, oscillations of the flux density. Same as in 
Figure [M 
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Fig. 60. — : Torsional mode, high magnetic field, oscihations of the flux density. Same as in 
Figure [191 
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Fig. 61. — : Torsional mode, high magnetic field, oscillations of the fiux density. Same as in 
Figure [201 
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Fig. 62. — : Torsional mode, high magnetic field, oscillations of polarization. Same as in Figure [211 
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Fig. 63. — : Torsional mode, high magnetic field, oscillations of polarization. Same as in Figure [191 
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Fig. 64. — : Torsional mode, high magnetic field, oscillations of polarization. Same as in Figure [20l 



